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We develop a three dimensional semiclassical theory which generalizes the Valet-Fert model in 
order to account for non-coUinear systems with magnetic texture, including e.g. domain walls or 
magnetic vortices. The theory allows for spin transverse to the magnetization to penetrate in- 
side the ferromagnet over a finite length and properly accounts for the Sharvin resistances. For 
ferromagnetic-normal-ferromagnetic multilayers where the current is injected in the plane of the 
layers (CIP), we predict the existence of a non zero mesoscopic GIF Giant Magneto- Resistance 
(GMR) at the diffusive level. This mesoscopic GIF-GMR, which adds to the usual ballistic contri- 
butions, has a non monotonic spatial variation and is reminiscent of conductance quantization in 
the layers. Furthermore, we study the spin transfer torque in spin valve nanopillars. We find that 
when the magnetization direction is non uniform inside the free layer, the spin torque changes very 
significantly and simple one-dimensional calculations cease to be reliable. 

PACS numbers: 72.25.Ba, 75.47.-m, 75.70.Cn, 85.75.-d 



Quantum effects in electronic transport are usually 
small deviations to the classical Ohm's law. Famous ex- 
amples include weak localization corrections^'^ and uni- 
versal conductance fluctuations^'^ which can be observed 
in diffusive systems. On the other hand, in small meso- 
scopic systems quantum mechanics can lead to an entirely 
different behavior of the conductance. Perhaps the most 
paradoxical example is the case of a smooth nanoconstric- 
tion in a two dimensional electron gas (quantum point 
contact) where by varying the strength of the confine- 
ment the conductance variation has a step-like character 
with plateaus quantized in unit of 2e? jlt''^ . This obser- 
vation was at first a bit puzzling since the quantum point 
contact does not have any source of scattering, neither 
elastic nor inelastic. The crucial concept for clarifying 
the picture was the notion of reservoirs (electrodes) at- 
tached to the quantum point contact where the energy 
relaxation takes place. 

We revisit bellow this issue in the context of the Gi- 
ant Magnetoresistance (GMR) effect^'* observed in mag- 
netic multilayers comprising ferromagnetic layers (F) sep- 
arated by a normal metal (N) structures. We consider 
F|N|F trilayer spin valve structure. There exist two ge- 
ometries for GMR. In the CPP geometry one injects 
the current perpendicular to the plane^^^^ of the lay- 
ers. The CPP-GMR can be understood simply within 
a two-current model where electrons with up and down 
spins experience different resistances as they cross the 
two magnetic layers. As a result, the configuration where 
the magnetizations of the FM layers are parallel (P) has 
a different resistance from the anti-parallel (AP) one, 
hence giving rise to the GMR. However, the original ex- 
periments^'®'^^'^'^ were performed in CIP geometry where 



one injects the current within the plane (CIP) of the lay- 
ers. The CIP setup is much simpler experimentally but 
the two current model (as well as its generalization, the 
Valet-Fert^"* drift diffusion theory) predicts a vanishing 
GMR. This usually complicates the interpretation of the 
experiments as one needs to introduce (sub mean free 
path quantum) microscopic approaches and the result- 
ing GMR can be quite sensitive to the details of the 
model (15^19 and references therein). Here, we predict 
the existence of an additional contribution to GIP-GMR 
which already exists at the drift-diffusion level. This 
contribution has the same origin as the quantification 
of conductance and should dominate the usual ballistic 
contributions for mesoscopic systems. 

In this work we address the following problems. First, 
a physical explanation for the role of conductance quan- 
tization in GIP-GMR is provided. Next, we develop a 
theoretical framework allowing for a quantitative predic- 
tion of this effect. This framework which we refer to 
as Continuous Random Matrix Theory in 3 Dimensions 
(GRMT3D) extends on a previous one dimensional semi- 
classical approach GRMT and goes beyond existing mod- 
cls^". Finally, we apply CRMT3D to the study of the spin 
transfer torque effect in a CPP nanopillar. 



I. SHARVIN RESISTANCE AND CIP-GMR: 
PHYSICAL PICTURE. 

A very transparent way to describe transport in a 
quantum system is the Landauer-Biittiker formalism^* 
where one describes quantum transport with a scatter- 
ing matrix that relates the amplitudes of outgoing to the 



2 



incoming modes. The quantum system is treated as a 
waveguide with A^^h ~ ^/^^ propagating modes. In the 
metalhc systems considered here, their cross section A is 
much larger than the square of the Fermi wave length 
Xp, so they typically contain thousands of conducting 
channels. The Landauer formula relates the transmis- 
sion probability r„ of channel n with the conductance as 
g = e'^ / hJ2n=i ^^^^ for a perfectly transparent sys- 
tem (T„ = 1) the system has a finite Sharvin resistance^^ 
3?sii = ^/(Ache^) . As mentioned above, the Sharvin re- 
sistance (also known as contact resistance) has been ob- 
served repeatedly, in particular in semiconductor based 
mesoscopic systems where A'ch can be tuned with the help 
of an electric field. The Sharvin resistance needs to be 
accounted for only once: when two systems A and B (of 
resistances "Ra and CRs) are connected in series, the total 
resistance is given (in the many channels limit A"ch ^ 1 
considered here) by^"^'^* "^ab = + — 3^sh- In other 
words one adds up the intrinsic resistances 3^^— 3^sh of the 
conductors in series (regular Ohm's law) and adds once 
a contact 3^sh/2 resistor for each electrode. For metallic 
spin valves under consideration here, the Sharvin resis- 
tance is the leading quantum correction to Ohm's law. 
3^sh typically corresponds to the resistance of an inter- 
face between two different metals or, say, 10 nm of bulk 
material. This resistance would normally be difficult to 
distinguish from a series resistance coming from the mea- 
suring apparatus. However, we shall see that in the case 
of the CIP geometry, the OMR signal simply vanishes at 
the purely classical level (more precisely in the limit of 
validity of the Valet-Fert model^^ described below) and 
the presence of the Sharvin (quantum) resistance pro- 
vides the leading source of GMR of mesoscopic samples. 
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Figure 1: Upper panels; Cartoon of the spin resolved chemical 
potential (up spins: dotted line, down spins: full line) of the 
top layer of a CIP-GMR setup as a function of the position 
X in between the electrodes. The chemical potentials drop 
from eV (x = 0) to (a; = L). The initial and final drop 
at a; = 0/L is due to the presence of the Sharvin resistances. 
Lower panels: numerical simulations of spin accumulation in 
a Co3|Agi|Co3 trilayer with a size L = 20nm. Left (a,c) and 
Right (b,d) panels correspond respectively to the parallel and 
anti parallel configurations. 



The proper generalization of Ohm's law to a three 
dimensional magnetic multilayer stack is given by the 
Valet-Fert equations: 



(1) 
(2) 



where Eq.(l) is Ohm's law relating the spin resolved 
current jo- to the gradient of the spin resolved chemi- 
cal potential Ha with the spin dependent resistivity pa. 
Eq. (2) expresses the (lack of) conservation of spin cur- 
rent: the divergence of spin current is balanced by spin 
relaxation which is proportional to the spin accumulation 
A/x = Hi and controlled by the spin diffusive length 

A cartoon of the system is presented in the lower panel 
of Fig. 1: it consists of two magnetic layers Fa and Fb 
separated by a normal spacer A^ and connected sideways 
to the two electrodes to which a voltage V is applied. To 
elucidate the role of Sharvin resistances in CIP-GMR, 
let us first ignore the role of spin-flip processes. In the 
absence of Sharvin resistances, one finds that the spin 
dependence of the resistivity in the magnetic layers is 
essentially irrelevant: the chemical potential must drop 
linearly from eV at x = to at x L irrespectively 
of the values of the resistivities pa- As a result, ix^ix, z) 
is constant along the growth z direction and there is no 
current fiow along z. There is neither spin accumula- 
tion in the system nor GMR. The same conclusion can 
be drawn from a full analysis of the Valet-Fert equations. 
The situation changes drastically when the system is con- 
nected in series with its contact (Sharvin) resistances as 
schematically sketched in Fig. 1 (a) and (b): for an elec- 
tron species (say majority electron) with low resistivity 
{paL <ti 3?sh) the chemical potential drops mostly at the 
contacts while for an electron species (say minority elec- 
tron) with large resistivity (paL 3> 3?sh) most of the drop 
takes place in the bulk. As a result some spin accumu- 
lation builds up in the system. In particular, in the AP 
configuration, the up spin (for example) chemical poten- 
tial varies along the z direction leading thereby to some 
spin current flow along the z axis. The current patterns 
become different for P and AP conflgurations and the 
GMR is restored. The color code of the lower panels 
of Fig. 1 represents the spin accumulation proflle calcu- 
lated A/x for a CoajAgilCoa (thicknesses in nm) trilayer 
using the theoretical approach described in the following 
section. We observe a clear non zero spin accumulation 
in Fig. 1 (c) and (d) for a CIP geometry. This effect 
is a direct consequence of the presence of non negligible 
quantum contact resistances. 



II. 3D SEMI-CLASSICAL THEORY : CRMT3D. 

In order to provide a quantitative description of afore- 
mentioned effects and to capture situations where the 



3 



(a) 




(6) 



© 



pout 



pout 



Figure 2: Panel (a) is a cartoon of discretization of tlie sys- 
tems in nodes connected by links. Each node is labeled by 
a latin index (i, j,...) and a link by the two latin indices of 
the associated nodes. Panel (b) is a cartoon of the Scattering 
matrix Sij. We define the region i on the left and j on the 
right. The in (out) stands for incoming (outgoing) probabili- 
ties mode Pi°^°"*, where the first index correspond to the side 
on which the probability current is defined. 



currents 



ij , mx 



(3) 



J 



are 4-vectors that encapsulate the current probabilities 
for majority (P^^ ^ ) and minority (P^^ ^ ) electrons as well 
as spin currents transverse to the magnetization of the 
layer (P^^ ^^^). The theory is defined by two fundamental 
equations relating the outgoing to the incoming currents: 
one for the links and one for the node. 

The link equation is identical to its counterpart in one 
dimension: 



pout 
■pout 



where the Scattering matrix Sij 



(4) 



(5) 



is composed of 4 x 4 transmission t, t' and reflection f, 
f' material dependent subblocks. The Sij matrix of a 
thin slice of material of width b is parametrized by two 
matrices A* and A*": 



magnetization has a non-trivial texture (domain walls, 
vortices...), we develop a full 3D semi-classical theory 
hereafter referred to as CRMT3D. This approach can be 
viewed as a generalization of the 3D Valet-Fert theory 
that properly accounts for Sharvin resistances and non- 
coUinear situations. In addition, CRMT3D can also be 
considered as a continuous version of the generalized cir- 
cuit theory^^ or equivalently of the random matrix the- 
ory developed in Ref. [26] . CRMT3D is a straightforward 
generalization of the recently developed CRMT (Contin- 
uous Random Matrix Theory) for unidimensional sys- 
tems^'*'^''. We refer to Refs. [24,27] for a full derivation 
of the one dimensional CRMT theory. For completeness, 
we recall below the basic objects of the theory before 
extending it to three dimensions. 

A schematic cartoon of the structure of CRMT3D is 
presented in Fig. 2(a) where the system is discretized into 
many nodes of small volume. Here we choose a simple 
Cartesian mesh but this choice is not compulsory. The 
nodes are connected to their neighbors by links. This set 
of nodes and links forms a circuit theory entirely equiva- 
lent to the so called generalized circuit theory^^'^^. The 
theory has four basic variables per link ij: P°"*, PJ", 
P°"' and P" where the labels in (out) refer to currents 
going from (to) the nodes while the index order ij (ji) 
indicates that the probability current is defined on node 
i (node j) side as sketched in Fig. 2(b). The probability 



t = 1 - A* 6, f = A^ 5. (6) 

Finally, the matrices A* and A'' of a ferromagnetic metal 
are parametrized by four independent parameters (F-]-, 
F;, Fsf, F„ix) and read. 



A* = 



A'' = 
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where d = 3 is the spatial dimension. As expected, for 
a unidimensional system (c? = 1) Eqs. (7,8) reduce to 
the one obtained in the CRMT approach, Eqs. (36,37) of 
Ref. [27]. The meaning of the parameters are also iden- 
tical and thus linked to 5 physical lengths : the spin 
resolved mean free path = l/Fo-, the spin diffusion 
length fgf = [4Fsf(F-|--|-F^)]~^/^; the transverse spin pen- 
etration length £± and the Larmor precession length ^l- 
The two latter are encoded into the complex number 
Fmx = i/£± + i/£'L- Note that although the role of and 
can lead to interesting physics, in the numerical simu- 
lations performed in this paper, we restrict ourself to sit- 
uations where these lengths are very small (sub nanome- 
ter) so that spin torques essentially develop at the normal 
metal-ferrromagnet interface. For a normal material, we 
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substitute in Eqs. (7), 

so that A* remains invariant upon arbitrary rotation of 
the spin quantization axis. A'' is obtained by fohowing 
the same procedure but with Tst replaced by — Fsf. The 
scattering matrices describing the interface between two 
metals are strictly identical to those developed in the one 
dimensional CRMT case to which we refer for their ex- 
pression (See section E of Ref. 27 for details). 

Although the natural variables of the theory are the 
probabilities Pj"''°"* , they are intrinsically related to the 
spin resolved chemical potential 



and the spin resolved currents 



(9) 



(10) 



flowing from i to j. The equation at the node is obtained 
by enforcing two conditions. First the chemical potential 
Hi depends on i only (and not on the link ij): 



(11) 



where Zi is the set of neighbors of node i. Second, the 
current is conserved at each node, 



•src 
Ji 1 



(12) 



where the source term ']f^ is present only for the nodes 
connected to electrodes. It is simply given by 



-IXT with, ixf 




(13) 



where Va is the voltage imposed at the electrode a. 
Rewriting Eqs. (11,12) in terms of probabilities and sub- 
stituting /ii, Eq. (9) yields the node equation, 
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TTTll^'^ (14) 



where \Zi \ is the total coordination number of the node i 
(counting the connections to other nodes plus the possi- 
ble presence of a connected electrode). The set of Eq.(4), 
Eq.(13) and Eq.(14) fully defines the theory. It is equiv- 
alent by construction to CRMT for one dimensional case 
and one easily verify that taking the continuous limit 
of Eqs. (9,10) for collinear system, one recovers the VF 
equations Eqs. (1,2). CRMT3D can be used in a variety 
of ways, both analytical and numerical. A very efficient 
numerical solution (used in the next section for up to a 
million nodes) consists of simply iterating the set of equa- 
tions (4). (13) and (14) from an arbitrary starting point 
until until convergence. 



III. NUMERICAL RESULTS FOR CIP-GMR. 

We now apply CRMT3D to CIP-GMR. Wc perform 
our simulations on various stacks on square samples of 
size L X L. A typical result is presented in Fig. 3 where 
the GMR defined as GMR= (31ap - :Rp)/{^ap) is plot- 
ted as a function of L for several F|N|F structures. Here 
Rp and Rap represent the resistance in the P and AP 
configurations, respectively. The GMR vanishes in two 
limiting cases: (i) when L <^ min(£-|-, £i), the resistance is 
entirely dominated by the Sharvin resistance which does 
not depend on spin; (ii) when L ^ max(^-|-, £j^) the re- 
sistance is dominated by the Ohmic resistance and spin 
accumulation vanishes as discussed above. Hence, one 
observes a negative correction for GMR (typically — 1%) 
for sizes L £^ (i.e. when intrinsic and Sharvin resis- 
tances have comparable contributions) . The actual value 
of the GMR depends on the kind of material considered 
(as shown Fig. 3) and the various thicknesses of the lay- 
ers. For instance, weakly resistive normal metals such as 
the copper Cu (blue squares and green triangles in Fig. 3) 
favors a shunting effect through the spacer which reduces 
the mesoscopic GMR signal. 




Figure 3: Numerical simulations of CIP-GMR (percent) as 
a function of the system size L for various F|N|F trilay- 
ers: Co3|Agi|Co3 (Black squares), Py3|Agi|Py3 (red circles), 
Py3|Cui|Py3 (blue diamonds) and Co3|Cui|Co3 (green trian- 
gles). 



In addition, we note the following characteristics: (i) 
We expect that raising the temperature has opposite ef- 
fects on the two sides of the negative peak. Indeed, when 
raising the temperature the effective mean free paths 
(which includes both elastic and inelastic scattering) de- 
creases. This makes the Sharvin contribution even less 
significant for large systems so that the GMR decreases 
in magnitude. However, for very small systems where 
the Sharvin contribution dominates, the GMR will start 
to build up. (ii) The sign of the effect effect is opposite 
to the CPP case: in the AP configuration, the minority 



5 



electrons of one layer take advantage of the CIP config- 
uration to propagate more freely in the other layer, (iii) 
One should keep in mind that this effect occurs in ad- 
dition to other microscopic ballistic contributions (typ- 
ically a few to ten %). Parametrically, this mesoscopic 
CIP-GMR vanishes algebrically as 1/(1 + L/ta) while 
ballistic contributions decay exponentially e'^^l^" with 
the width of the normal spacer. It is therefore possi- 
ble to see the mesoscopic effect only but for most stacks 
one would observe both effects simultaneously and meso- 
scopic CIP-GMR should therefore be observed as a dip 
in the L dependence of GMR. 

Measuring the size dependence of CIP-GMR is not an 
easy experimental task. However, a very similar signa- 
ture can be obtained by measuring the (two-terminal) 
GMR using a STM tip as a function of the distance be- 
tween the tip and the contact electrode. The setup is 
presented in the inset of Fig. 4 where the tip is placed on 
top of the stack at a voltage V while the two other elec- 
trodes are grounded. It corresponds to a geometry and 
sizes very similar to those used in the low temperature 
STM experiment of Ref.^*. As the distance x between 
the tip and the contact electrode increases, the GMR is 
anticipated to change from positive (GPP like) to nega- 
tive when the mesoscopic CIP GMR effects dominates. 
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Figure 5: Torque r per total current / as a function of the 
local angle 9 between the magnetization of the free layer 
and the reference magnetization of the thick layer for a 
Cu48|IrMn8|Py8|Cu4|Py8|Cu2o|Ptg|Cuioo pillar for a current 
density I = 2.5 A.cm~^ . The full black line corresponds to 
ID CRMT calculation while the cloud of red dots corresponds 
to the CRMT3D calculation. Inset: reference magnetic con- 
figuration of the free layer, calculated with micro-magnetic 
simulations^'' . 
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Figure 4: mesoscopic GMR in percent as function of the STM 
position X for a F|N|F trilayers: Co2|Ag2|Co2, with a length 
L = 500nm and a width W = 60nm. Inset: Cartoon of the 
top view of the setup. 



IV. CONCLUSION: SPIN-TORQUE IN A 
REALISTIC CPP SPIN VALVE. 

To conclude this paper, we take advantage of the ca- 
pability of CRMT3D to treat systems with magnetic 
texture and perform a study of spin transfer torque in 



a spin valve. We aim at evaluating the role of mag- 
netic texture which is often disregarded in the calcu- 
lation of spin torque made in micro-magnetic simula- 
tions. Our nanopillar corresponds to the following stack: 
Cu48|IrMn8|Py8|Cu4|Py8|Cu2o|Pt8|Cuioo- It consists of 
a polarizing layer pinned by the IrMn antiferromagnetic 
layer and a free 8 nm permalloy layer. This setup basi- 
cally corresponds to the experiments reported in Rcf."^^ 
and has been designed in such a way that the current in- 
duced magnetization reversal behaves in a coherent way 
(i.e. as close to a macrospin as possible) so that in this sit- 
uation the role of magnetic texture is believed to be fairly 
small. Nevertheless, the Oersted field which is present at 
high current introduces a small "banana shape" magnetic 
texture, as shown in the inset of Fig. 5. Our starting 
point is the corresponding stationary magnetic configu- 
ration of the free layer obtained from a micromagnetic 
simulation in presence of the Oesterd field^^ . In a second 
step, we perform two different calculations of the spin 
torque: (I) a full CRMT3D calculation the local spin 
transfer torque ti{x, y) in presence of the Banana shape 
magnetic texture (the polarizing layer which is pinned by 
the IrMn layer is supposed to have no magnetic texture) . 
(II) Wc take an approach which ignores the role of in 
plane spin currents: one assumes that current density is 
homogeneous across the nanopillar and parametrize the 
local spin transfer torque r = /(6') as a function of the 
angle 9 between the (local) magnetization and the ref- 
erence fixed polarizing layer. When the system acquires 
some magnetic texture one uses Tij{x,y) = f{9(x,y)). 
The paramctrization f{9) is obtained using the one di- 
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mensional version of CRMT. The effective ID approach 
(II) is equivalent to the 3D approach (I) in the absence of 
magnetic texture and has become quite common in dy- 
namical micromagnetic studies of current induced phe- 
nomena'^^"^^. The results are shown in Fig. 5 where the 
effective ID CRMT approach t = f{9) (line) is con- 
trasted with the full 3D calculation where Tii{x,y) is 
plotted as a function of 9{x, y) (red dots). The apparent 
" noise" of the CRMT3D calculation reflects the fact that 
the torque is not a function of 9 only but fully depends on 
the spatial position {x,y). One can see that even though 
the general picture is captured by the effective one dimen- 
sional approach, a typical error of more than 10% may be 
observed. We expect that upon performing an integra- 
tion of the (highly nonlinear) 3D Landau-Lifshitz-Gilbert 
equation, such a systematic error will result in strong in- 
accuracy, even in the favorable situation considered here. 



We conclude that micromagnetic simulations of real pre- 
dictive power, which are highly desirable for spintronic 
applications, require to treat magnetic and transport de- 
grees of freedom on equal footing. In particular, a natural 
route would be to perform full CRMT3D calculations of 
the spin transport properties of the device "on the fly" 
during the micromagnetic simulation. 
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